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An influence -function technique is employed to analyze the stresses 
and deformations due to flaws in an idealized composite material com- 
posed of parallel, equally spaced, tension -carrying filaments embedded 
m a shear -carrying matrix. Static tensile stress concentration factors 
are obtained for the cases of two equal-length transverse collinear cuts 
across various numbers of filaments, and for the case of a periodic 
array of transverse collinear cuts. In addition, a static shear stress 
concentration factor is obtained for a single cut across an arbitrary 
number of consecutive filaments. Dynamic stress concentration factors 
are obtained for two cases in which two collinear cuts are suddenly 
introduced into a stretched filamentary sheet. Matrix shear loads are 
investigated and their variation in the longitudinal direction is studied 
for several cases involving single and double cuts. In addition, loads 
m broken filaments are calculated for some single-cut eases, and their 
implications for an existing statistical tensile failure analysis for 

i * v 

composite materials are briefly discussed. 

In the case of two collinear cuts it is found that interaction 
between the two cuts is significant only when the distance between the 
cuts is no greater than the cut length and, hence, that more widely 
spaced cuts may, for all practical purposes, be treated as isolated 
cuts. The interaction between closely spaced cuts, however, is pro- 
nounced. It is found that two closely spaced cuts can cause higher 
stress concentrations than a single cut across a comparable total 



number of filaments. This result suggests that a design criterion 
"based on the residual strength of a composite weakened by a single cut 
of prescribed length can be unconservative, even though only cuts of 
lesser length are present. 

The stress concentration factors for a periodic array of collmear 
cuts, which for the present model are equivalent to the stress concen- 
tration factors for a transverse cut in a strip of filament- stiffened 
material, are found to agree closely with results of an earlier approxi- 
mate analysis, except for cut lengths approaching the width of the 
strip. This indicates that m the practical range of interest of the 
cut-length to strip-width ratio, the approximate stress concentration 
factor is sufficiently accurate for engineering purposes. 

The matrix shear stress concentration factor for a single cut of 
arbitrary length is seen to increase with cut length somewhat more 
rapidly than the tensile stress concentration factor, particularly for 
the shorter cuts. This suggests that some filamentary composites 
might be primarily susceptible to tensile failure when weakened by 
small flaws, hut might be more susceptible to matrix shear failure 
when larger flaws are present. 

Results of limited calculations of filament dynamic stress con- 
centration factors for two collmear cuts are m agreement with earlier 
results for a single cut, and support the previous conclusion that 
dynamic effects are of secondary importance for the type of model 
investigated. 



Calculations of the longitudinal variation .of matrix shear forces 
show that they decay less rapidly with increasing cut length. Broken 
filaments are seen to recover load less rapidly with greater cut length, 
a result which suggests that tensile failure analyses for imperfect 
composites, depending strongly on the ineffective length of broken 
filaments, should account for the variation of ineffective length with 
flaw size, rather than employ a single value of ineffective length 
regardless of the distribution of initial imperfections. 
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INTRODUCTION 


Composite materials are finding increasingly wide application in 
aerospace structures. For example, many solid-propellant rocket motor 
cases are constructed "by winding resin-coated glass filaments on a 
mandrel. Coated fabrics, because of their great flexibility, have been 
used m applications requiring the temporary packaging of large, low- 
density structures into small storage volumes Currently, much effort 
is being applied to the development of lightweight composite materials, 
which typically are composed of high-modulus filaments embedded m 
plastics or metals of relatively low density and modulus. 

The rational design of a structure requires knowledge of the 
stresses which it is likely to experience. However, because of their 
inhomogeneity, composite structures often do not lend themselves to 
representation by tractable mathematical models. For example, the 
walls of filament -wound rocket motor cases are constructed of numerous 
layers of windings, with the winding direction varying from one group 

\ t 

of layers to another. A structure of such complexity presents 
formidable analytical difficulties. However, if attention is focused, 
on a single layer of composite material, it is possible to obtain 
analytical results which may be applicable ultimately to more complex 
arrangement s . 

The problem of stress concentration m aerospace structures is one 
of continuing importance. In reference 4 , a stress concentration 
problem for a plane of parallel, equally spaced filaments embedded m 
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a matrix was formulated and solved by Hedgepeth. The principal results 
of the analysis were the static and dynamic stress concentration factors 
caused by a single cut across a number of adjacent filaments. In 
reference this analysis was extended by Hedgepeth and Van Dyke to 
problems involving some special two-dimensional distributions of 
parallel filaments. 

In this thesis, the analysis of reference k is extended to the 
cases of two coHmear cuts and of periodic collmear cuts of equal 
length. Also, the single-cut problem of reference 4 is re-examined 

T T 

for the purpose of more fully exploiting the potential of the model * 1 

4 

for the study of some composite materials. Specifically, matrix shear 

» 

forces, and the decay of matrix shear forces and filament tensile 

forces with axial distance from the cut are studied. In addition, 

* . 

an expression for maximum shear force, analogous to the filament stress 
concentration factor, is determined for the case of a single cut. 

In the case of a double cut, static stress concentration factors 
are computed and compared with some related single -cut results. In 
addition, dynamic stress concentration factors are calculated for two 
of the simplest double-cut cases. 

For the problem of periodic collmear cuts, which is related to 
the problem of a finite -width sheet weakened by a central cut, static 
stress concentration factors are calculated for various combinations 
of cut length and distance between cuts, and are compared with results 
of an approximate analysis for an isotropic elastic sheet 
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REVIEW OF LITERATURE 

An extensive survey of research m composite materials, with 
particular emphasis on fibrous or filamentary composites, is contained 
in reference 6, where numerous additional references are cited. 

Particularly relevant to the present work are references 4 and 5- 

% 

In reference 4, an analysis is presented of the stress con- 
centration around a single straight cut across an arbitrary number of 
filaments in an idealized composite composed of a single infinite 
layer of parallel, equally spaced, tension-carrying filaments embedded 
in a shear-carrying matrix. A closed-form expression was obtained for 
the stress concentration factor (defined as the ratio of the highest 
load in an unbroken filament to the far-field applied load) as a 
function of the number of broken filaments. This result was substan- 
tiated experimentally by Zender and Deaton m reference 9. In 1 ' 
reference 4, other phenomena of interest, such as matrix shear forces, 
and the variation of filament and matrix loads in the filament 
direction, were not investigated. 

In reference the influence function technique introduced m 
reference 4 was extended to the problem of stress concentration in a 
composite material composed of two-dimensional arrays of filaments 
embedded in a shear-carrying matrix. Because of the difficulty of 
obtaining closed-form inversions for various transformed quantities, 
extensive numerical computation was required. It was found that, in 
general, the breaking of a given number of neighboring filaments causes 



k 


greater stress concentration in the single layer than in the two- 

dimensional array. An additional analysis of a single layer , m which 

the matrix material around a single broken filament was assumed to 

undergo ideal plastic deformation over a portion of the filament 

length, indicated that the inclusion of matrix plasticity mitigates ‘ 

the filament stress concentration effect predicted by the ‘elastic 

analysis. However, it appears that the results are too limited for 

sweeping conclusions to be drawn regarding the effects of plastic 

» 

deformation m an imperfect filamentary composite. Experimental 

* 

verification of the results of reference in view*" of the necessarily 

•* 

more complicated nature of the mathematical models, is .likely to prove 
much more difficult than was the case with the single-layer elastic 


model 
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GOVERNING EQUATIONS 

I 

As far as is practical the notation employed in reference ij- is 
retained m the present analysis. The configuration is shown in 
figure 1, along with the coordinate system and some notation. The 

i 1 

analytical model is one which is commonly used in "shear-lag" analyses. 
It is composed of an infinite -single layer of parallel tension -carrying 

- . i 

members (filament) embedded in a matrix which carries only shear. The 
filaments are separated by a constant distance d and are numbered 

t 

from -oo to to from the bottom upward. The coordinate along the 
filaments is x and the displacement of the n"^ filament at location 
x is u n (x, t) . The force m the filament is denoted by p Q (x,t) 

and is given m terms of u n by 




where EA is the extensional stiffness of the filament. The shear 
force per unit length between the n^ 1 and (n+l)^ 1 filaments is 
defined here by s n = ^(u n+1 - u n ) . Conservation of momentum of the 
n th filament then requires 




where the assumption has been made that the mass per unit length j 
associated with the n^ 1 filament is concentrated at that filament. In 
terms of displacements, the equation of motion becomes 




( 2 ) 
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M % + ' 2 "» + ^ ‘ r % 


In figure 1, filaments -2, -1, 2 and 3 are shown broken at x = 0, 
with the remaining filaments intact. In general, for two collinear cuts 
(at x = 0) through q and s filaments, let - (m + q) ^ n % - (m + l) 
and r + l^n = r + s, respectively, denote the broken filaments, the 
two cuts being separated by r + m + 1 intact filaments. Then the 
boundary conditions are 

u n (0,t) =0, - m ^ n s r, - (m + q + l) = n, n = r + s + 1 

(5) 


and 


p n (0,t) = 0, - (m + q) 5 n 5 - (m + l), r+l<n5r+s 

For x large, the force in each filament approaches the uniform applied 
force, denoted by p; that is 

p n (±«>,t)=p (*0 


For the time -dependent problem, the following initial conditions are 
required: 

p n (x,0) = p 


du n 

(x,0) = 0 

Bt 


(5) 
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For convenience , let 

Pn = P p n 

^ ” P ^EAGh Un 
s = p (JO: Sq 

VEAd ^ ( 6) 

x = 

V Gh 

t = IEr 

~ \j Gh T 

Then the equilibrium equation becomes 

8%„ o 2 b„ 

— - + u n+1 - aj„ + i^-i r (7) 

d | 2 dr 2 

with boundary conditions 

U n (0, T ) =0, ■ i 5 n ^ r, - (m + q + 1) ^ n, n^r + s + 1 
P n (0,r) = 0, - (m + q) = n = - (m + l), r + 1 ^ n 5 r + s (8) 
P n (*»,T) = 1 



and initial conditions 
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SOLUTIONS 


Stress Concentration Factors 

The ‘boundary value problem for static loading is defined by 
equation ( 7 ) with the right hand side set equal to zero, and boundary 
conditions ( 8 ). The solution is complicated by the fact that the 
boundary conditions at £ = 0 are mixed. This difficulty was overcome 
m reference 4 by use of an influence function technique. For a 
thorough discussion of this technique, the reader is referred to 
reference 4. The influence functions V n (5) and N n ( 5 ) = dV n (|)/d| 
are, respectively, the nondimensional displacement and force in the 
nth filament when the filament sheet is completely cut along 5=0 
and the zeroth filament is displaced axially a unit amount at 5 = 0 
while all other filaments are held fixed at 5=0. In terms of V n (5) 
and N n (5), the dimensionless force and displacement are given by 


00 

p nU) = 1 + £ N n .i(6)Ui(0) 

i=-co 

CO 

U n (5) = § + Y V n _ i (5)U i (0) 

l=-eo 

In reference K V n ({) was found to be 


(n) 


V n (5) - 



cos n 0 e 


■ 2 |sin 0/2 


d 0 


( 12 ) 


Application of boundary conditions ( 8 ) to equations (ll)- yields first 
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-(m+1) 


r+s 


P n (£) = 1 + Y, n n-i^) u i(°) + Y K n-i^) u i(°) 

i=- (m+q) 


-to) 


i=r+l 

r+s 


(13) 


u nU) = I + Y V n-i(^ u x(°) + Y Vi(5) u i(°) 

i=- (m+q) 


i=r+l 


since U^O) =0 for other values of i, and second 

-(» + 4 )SnS-(m + X) 


0 = 


H v~> -(.m + q; ^n ^ ■ 

1 + ) W n -i(°) u i(°) + ^ w n-i(°) u i( 0 )^ and 

i=-lfm+q) i=r+l r + l'5n 5r + 


to 


■which expresses the condition of zero load on the ends of the broken 
filaments. 

Equations (l4) constitute a set of q + s linear algebraic 

equations in the q + s unknowns, U n (0). Their solution set can be 

substituted into equations (13) to obtain expressions for load and 

* 

displacement m any filament. However, before this can be done, the 

‘ 

integral representation of the influence functions", V n (f ) 'and H n (i) 
must be evaluated. Only H n (0 ) = dV£(0)/di was .evaluated m 

reference k because the computation of stress concentration factors, 

\ ' 

which was the main purpose of that investigation, does not require the 
evaluation of the influence functions for non-zero values of £. 
Equation (12) can also be written as 

V n U) =~ J cos 2n0e“ 2 ^ Sin0 de 


(15) 



Ik 


Integration of this expression has been carried out in Appendix A. The 
result is 


V n (£) = (-1) 


n 


n-1 


l2n(2£) - L2n(25) - 1 £ (-1) 3 

k=0 


2n-2k-l 


* r (*' * £) 

r(sn + i - k) 


I# n = 0 


(16) 

in 'which I p n is the modified Bessel function of the first kind, L p n 
is the modified Struve function (see reference 1), and the finite sum 
must be taken to be zero for n = 0. A similar expression for V n (i) 
can be derived for n < 0; this is not necessary, however, since it can 
be seen from equation (15) that V_ n (|) = V n (|). Differentiation of 
equation (16) yields 


N n (|) = (-l) n 


+ X 2n-l( 2 ^) - < L 2n+l(2t) + ^^(Sg) 





(2n - 2k - l)r^k + 

r ( 2n + I " k ) 


2n-2k-2 


n 


£ 0 


( 17 ) 


in which the finite sum again is to be omitted for n = 0. From 
equation (17) it is found that 
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N n (o) 


i)- 

3t(4n^ - 1) 


(18) 


m agreement with reference 4. 

With substitution of the appropriate values from equation (l8), 
equations (l4), which are merely linear algebraic equations for the 
displacements of the ends of the broken filaments, can be solved. 

Then the solutions to equations (1*0, along with the influence functions 
given by equations (l6) and (17), can be substituted into equations (13) 
to obtain the load and displacement of any filament in the sheet. 

Tensile stress concentration factors for two collinear cuts .- In 
the case of a single cut in a filamentary sheet, the stress concentra- 
tion factor depends only upon the number of consecutive filaments 
traversed by the cut. In reference 4, the stress concentration factor 
for a cut across n consecutive filaments was found to be 


- 4*6»8» • • (2n + 2) 
“ 3*5*7* • • (2n + 1)' 


1,2,3, 


(19) 


where K n is the ratio of the maximum load in either of the two intact 
filaments directly adjacent to the cut, to the load at infinity. 

In the case of two cuts, two additional parameters appear: (l) the 
number of broken filaments in the second cut, and (2) the number of 
intact filaments separating the two cuts. In general, each cut may 
sever any number of filaments and any number of intact filaments may 
separate the two cuts. Hence the general case of two collinear cuts is 
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without symmetry, and calculations covering reasonably wide variations 
of the three pertinent parameters would entail a considerable computa- 
tional effort. However, if the problem is simplified by requiring the 
cuts to be of equal length, then the most essential features of the 
two-cut problem are retained 'while the computational effort is greatly 
reduced through consideration of the resulting symmetries. 

In what follows, then, it is assumed that the two cuts traverse 
the same number of filaments. In the analysis of this reduced problem 
one of two cases arises, depending on whether the number of intact 
filaments between the cuts (henceforth called "interior filaments") is 
even or odd. The analysis of these two cases is presented in 
Appendix B. 

The double-cut stress concentration factor K n ^ m is defined as the 
ratio of the greatest load m the most highly stressed unbroken filament 
to the load at infinity, for two collinear cuts, each across n fila- 
ments, and separated by m interior filaments. For an odd number 
(2r + l) of interior filaments, the stress concentration factor for 
two collinear cuts, each across n filaments is found in Appendix B 
to be 

r+n 

^n, 2r+l = 1 + Yj 
i-r+1 

where each 1^(0) is calculated from equations (b 4) . For an even 
number (2r) of interior filaments. 


r+i 


(°) + Vjo) 


u i(°) 


(26) 
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r+n-1 , ^ 

Kn,2r = 1 + L (W°> + H r .i. x (0)J U x (0) 


( 21 ) 


where each U^(0) is calculated from equations (BIO) . , * 

The stress concentration factors, given by equations (20) and (21), 

V ' k 

1 

have been computed for two equal-length collinear cuts, each traversing' 

* i w 

t 

from one to eight filaments and separated by intact filaments ranging in 
number from one to 1 6 . The results are presented in Table 1, and are 

plotted in figure 2 for the various values of n, the- number ‘of fila- 

( 4 < 

% 

ments severed by each cut. Although the curves are meaningful only for 

J t 

* J i * 

integral values of n, they are plotted, as, continuous* curves for illus- 

} » * 

^ t ( * 

trative purposes. Also shown in figure 2 >are the single-cut stress 
concentration factors K n , which are the asymptotic values of K n ^ m 
for large m. 

As can be seen in figure 2, the interaction between cuts is 
essentially a local phenomenon, being confined to separation distances 
(values of m) on the order of the cut length. The results indicate 
that cuts -which are separated by distances greater than their length 
may, for all practical purposes, be treated as isolated cuts. 

For closely spaced cuts (small m), however, the interaction is 
pronounced. Hence, the values of K n ^ m for small m are of particular 
interest, since they are associated with states of high stress concen- 
tration. In order to present a comparison of the severities of single 
and double cuts, the stress concentration factors for two cuts of 
length n separated by one and two intact filaments, and ^ 
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Table 1. Filament stress concentration factors 
for two collinear cuts of equal 
length 


Wo. of 
filaments 
between 
cuts, m 


Wo. of filaments 

in each cut, n 


l 

2 

3 

4 

5 

6 

7 

8 

1 

1.714- 

2.359 

2.964 

3.543 

4.102 

4.646 

5-177 

5.698 

2 

1.4-12 

1.788 

2.142 

2.481 

2.808 

3.126 

3.436 

3.740 

3 

1.368 

1.690 

1.989 

2.272 

2.543 

2.805 

3.061 

3-310 

4 

1.353 

1.654 

1.928 

2.185 

2.430 

2.666 

2.895 

3.118 

5 

1.346 

I.636 

1.897 

2.l40 

2.370 

2.590 

2.803 

3.009 

6 

1.342 

1.626 

1.879 

2.113 

2.333 

2.542 

2.745 

2.940 

7 

1.340 

1.620 

1.867 

2.095 

2.308 

2.510 

2.705 

2.893 

8 1 

1.338 

1.615 

1.859 

2.082 

2.290 

2.488 

2.677 

2.859 

9 

j 

1.612 

1.854 

2.073 

2.278 

2.471 

2.655 

2.833 

10 



1.849 

2.066 

2.268 

2.4-58 

2.639 

2.813 

11 



1.846 

2 06l 

2.261 

2.448 

2.626 

2.798 

12 



l 

2.057 

2.255 

2.44o 

2.616 

2.785 

13 




2.054 

2.250 

2.434 

2.608 

2.775 

l4 





2.246 

2.428 

2.601 

2 766 

15 





2.243 

2.424 

2.595 

2-759 

16 






2.420 

2.590 

2.753 
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Figure 2. Filament stress concentration factors for two equal— length 
collineax cuts 


are compared with K pn , K/^n+i* and K 2n +2 dn f^^nre 3* These parti- 

cular values of the single-cut stress concentration factor have been 
chosen because they represent either the , same total number (2n) of 
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broken filaments as does K n ^ m , or the number (2n +1 or 2n + 2) 
which would be broken if the pair of cuts separated by one or two intact 
filaments, characterized by or K n ^ 2 , were to coalesce. In 

figure 3 it can be seen that the greatest stress concentration factor 
is not necessarily associated with the greatest total number of broken 
filaments. This is especially true in the case of the longer cuts, 
due to the fact that the stress concentration factors for double cuts 
increase more rapidly with cut length than do the factors for single 
cuts across a comparable number of filaments. The results suggest that 
a design criterion based on the residual strength of a structural 
component weakened by a single cut of prescribed length can be uncon- 
servative, even though only cuts of lesser length are present. 

Shear force concentration factor for a single cut .- A result which 
was not presented in reference 4, hut which can he extracted from the 
analysis is the magnitude of the most severe matrix shear force due to 
a single cut of arbitrary length. In the case of a single cut the most 
severe matrix shear force occurs in the neighborhood of each end of the 
cut. From reference 4, for a cut which starts at the zeroeth filament 
and severs n filaments, the load and displacement of the 1 th 
filament are given hy 




2 3 4 5 6 7 8 

Length of each double cut, n 


Figure 3* Comparison of filament stress concentration factors for single 
and double cuts 
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p ± a) 


i + 


n-l 

I 

r=0 


N i-r (l)U r (°) 


( 22 ) 


n-l 

U ± (5) = £ + Y, v i-r^) u r(°) 

r=0 


(23) 


with 


n-l 

0 = 1 + X %_ r (0)U r (0), 0 = l = n - 1 (24) 

r=0 

The most severe (nond tmensional ) shear force is given by S n -l(°) (see 
equation (23) and the second of equations (10) ), -which is merely 


since U n (0) = 0 because the filament is intact, and U n _i(0) is 

the nondimensional displacement of the end of, the (n-l) filament. 

* 

From symmetry considerations .this can be written as 

s A-l(°) “ - u o(°) 

■ } * 

or 

W ^ | u o(°)| ' 

where Sj,^ is defined as the peak magnitude of the most severe matrix 
shear force. Of course, U o (0) varies with the number of broken 
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filaments and is found by solution of equations (24) for the appropriate 
value of n. 

Solution of equations (24) for the first six values of n yields 
the following values of S^y. 


n 

^ ^max 

1 

i 

1 

2 

3/2 

5 

15/8 

4 

35/16 


315/128 

6 

693/256 


Inspection of the first six values indicates that they conform to the 
expression 


jSmax 1.2-3.— <26) 

' “ 2 2n-2[(n - l)'] 2 

This formula has been shown in Appendix C to hold for all positive 
values of n. The quantity ^ is seen to be the ratio of the 

maximum shear force for n > 1 to the maximum, shear force for a 
single broken filament. Hence, it can be viewed as a matrix shear 
force " concentration factor," m the sense that it describes 
the growth of the maximum shear force with increasing cut length. 

In this context, its comparison with the filament stress concentration 
factor for a single cut (see equation (19) ) is appropriate. Therefore, 
these two quantities have been plotted in figure 4 for values of n 







up to 13. Both factors are unbounded as n approaches infinity. 
However, their relative magnitudes for large n can be determined by 
■writing equation (19) in the form 
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K n 


2 gn n! (n + 1) » 
(2n + l)t 


and forming the ratio 

^ ^a-x _ n(2n + l)[(2n)Q 2 
K n 2 2to " 1 (n + l)(n') 4 


(27) 


By the use of asymptotic formulas for the factorial functions, it is 
found that 


-> 1, as n -» co (28) 

% 

Equations (27) and (28) indicate that the dimensionless maximum matrix 
shear force and tensile force in the filaments are of comparable 
magnitude for all values of n; however, as seen in figure 4 , the shear 
force initially increases more rapidly with cut length than does the 
filament tensile force. This result suggests that some composites 
might be more prone to tensile failure when weakened by small flaws, 
but more susceptible to matrix shear failure when larger flaws are 
present. 

Tensile stress concentration factors for periodic coll inear cuts .- 
The periodic collinear cut configuration is illustrated in figure 
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Figure 5. Coordinate and notation systems for periodic collinear cuts 
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Each cut is across n filaments. The period of the cuts, that is, the 
distance (in number of filaments) between the first broken filaments in 
succeeding cuts, is denoted by the integer R. The periodic configura- 
tion is of interest because, within the limitations of the present 
model, the periodic -cut problem is equivalent to the problem of a 
central transverse cut across n filaments m a filamentary strip which 
is R filaments wide. Of course, for the problem to be meaningful, 

R must be greater than n. The expressions for load and displacement 
are given by equations (ll) subject to the mixed boundary conditions 

P-lCo) = 0, i = Rr, Rr + 1, ..., Rr+n-1, r = 0, ±1, ±2, ... 

(29) 


and 


U^(0) =0 for all other l 


where, for convenience, one of the cuts is begun at the zeroeth fila- 
ment. Application of boundary conditions (29) to equations (ll) yields 


P ± (S) = 1 + 


E ■- 

m=-«> v 


kJiKWo) +N i _ Km _i(l)u Rm+1 (o) 


+ 


+ %-Rm-n+l 


( £ ) u Rm+n-l 



.(50) 



20 


Uid) = I + £ {vi.M(S) D an(°) +T i-Km-l^) u an + l(°) 

m=-oo ^ 


+ . . . + V^_-g m _ n+ ^ ( £ )UR m 4- n . 


(3D 


since U^(o) = 0 for other values of i, and 

00 r 

0 = 1+ ^ | w i-Em(0) u Em(0) + %-Rm-l(°) u Rm+l(°) 
n=-a> '■ 

+ ••• + N i-Sm-n+l(°) u Rm+n-l(°)^ 7 Br = i = Br+n-l, 

r = 0, ±1, ±2, . . . 


(32) 


which expresses the condition of zero load on the ends of the broken 
filaments. Replacing l by Rr + q in equations-* (32) yields 


0 = 1 + 


00 f 

]T [%(r-m)+j(°) u Bm(°) + %(r-m)+j-l(°) u Rm+l(°) 

m=~oo v 


+ ... + % ( r -m) +j -n+1 (° ) u Em+n-l 0 ~ 3 ,-h " 1 


Now letting k = r - m gives 


0 = 1 + ^ I ^Rk+j (° ^R(r-k) ^ + n Rk+a-l^°) U R(r-k)+l(°^ 
k=-oo v- 


• • • 


+ N Rk+j-n+l(°) u R(r-k)+n-l(°) f > 0 = j = n - 1 (33) 
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However, due to the periodicity of the displacements, 

u R(r~k)(°) = u o(°)-> u R(r-k)+l = U x (0), etc., so that equations (33) 

become 


° = 1+ Z { H Rk+ 0 (0)U o(°) +w Rk+j-l(°)Ul(0) 

k=-a> 

+ ••• + N Rk 4 -j-n+l(°)Vl(°)J > 0 = 3 = n - 1 

which can also he written as 

00 00 

° = l+U o (0) Z H Ek+J (0) +U 1 ( °) Z N Rk + j-i (0) 

k=-co k=-oo 

00 

+ ... +U n _ 1 (0) £ NRk+ 5 -n+l(0), 0 - 0 - n - 1 (34) 

k=-co 

which are merely n linear algebraic equations for the (0) ' s. 

It remains to evaluate the coefficients given by the infinite series. 
This has been done in Appendix D, and the results are 



A further slight simplification of equations (30), (31)> and (34) is 
afforded by noting the symmetry of the displacements ' m each cut, which 
gives U Q (0) U n _;j_(0), Uj_(0) = U n _2(0), etc. This simplification 
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reduces the order of the system lay roughly half (actually, by n/2 

n-1 , 

equations for n even, and by — for n odd). 

The filament stress concentration factor is obtained by calculating 

the maximum load m the first intact filament at either end of any of 

the cuts. For this purpose, the n^* 1 filament is chosen; thus the 

stress concentration factor is given by 

n even. 


co .. 

Kn = p n(°) - 1 +U o (0) ^ |ftn-Rk(0) + Ni_Rk(0)J> 


Vl(°) £ {>VBk-l( 0 ) + %-Kk(0)j + ••• 
U n/2 -1 (0) I fn/2 W°> + \/S -Ek<°> 

k=-co 


n odd: 
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where it should be noted that a* 1 *! 'where the 

U 1 (0) , s have been calculated from equations (3*0 for specified values 
of n and R (cut length and period). 

The stress concentration factor has been calculated for 

values of n from 1 to 6, and for values of R up to 36. The results 
are presented in Table 2. The results are plotted in figure 6, for a 
limited range of R - n, in terns of net section load (average load in 
the unbroken filaments) rather than average load at infinity, in 
keeping with customary engineering practice. Again, although the 
results are meaningful only for integral values of R - n, they are 
plotted as continuous curves for ease of illustration. Also plotted in 
figure 6 are results of an approximate analysis by Dixon (see refer- 
ence 3); which were employed by Zender and Deaton in reference 9 to 
convert the infinite- sheet results of Hedgepeth (reference 4) to a 
form usable for analyzing their experimental data on strips of fila- 
mentary material. As can he seen in figure 6, the results of the 
approximate analysis of Dixon are in close agreement with the present 
results except for cuts whose length approaches the width of the strip. 
In this range, the predictions of reference 3 are seen to underestimate 
the stress concentration factor. However, in the range of cut length 
to strip width ratios of practical interest, Dixon’s approximate 
formula, although slightly unconservative, is quite sufficiently 
accurate to justify its use. For example, although the comparison is 

not shown here, Dixon' s results, when applied to the data of Zender 

» 



Table 2 . Filament stress concentration factors for 
periodic collmear cuts (based on net 
section load) 


Het 

Cut length, n 

DCtli .LULL y 

R-n 

1 

2 

3 

4 

5 

6 

1 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

2 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

3 

1.061 

1.085 

1.098 

1.105 

X • 110 

1 » X 13 

4 

1.106 

1.155 

1.182 

1.199 

1.210 

1.218 

5 

1.138 

1.209 

1.250 

1.277 

1-295 

1.308 

6 

1.163 

1.251 

1.305 

1.3^2 

1.367 

1.386 

7 

1.182 

1.285 

1.351 

1.396 

1.429 

1-453 

8 

1.198 

1.314 

1 390 

1.443 

1.482 

1.512 

9 

1.210 

1.337 

1.422 

1.483 

1.529 

1.563 

10 

1.221 

1.357 

1.451 

1.518 

1.569 

1.609 

11 

1.229 

1.37^ 

1.475 

1.549 

1.606 

1.650 

12 

1.237 

1-389 

1.497 

1.576 

' 1.638 

I.687 

15 

1.254 

1.424 

1.548 

1.642 

1.717 

I.778 

20 

1.272 

1.462 

1.605 

1.718 

1.810 

H- 

1.886 

25 

1.284 

1.487 

1.643 

1.769 

1.873 

1.961 • 

30 

1.291 

1.504 

1.670 

' 1.806 

1.919 

2.016 
















and Deaton (reference 9 )> were found to be indistinguishable from the 
present results. 

Dynamic stress concentration factors for two collinear cuts .- In 
reference k, dynamic stress concentration factors were found for cases 
in which a cut is suddenly introduced in a stretched filament sheet 
and, in a separate analysis, an apparent upper limit of 1.27 was found 
for the dynamic response factor (ratio of maximum dynamic to static 
stress concentration factor). Investigated were cuts across one, two, 
and three filaments, and the limiting case of a finite-length slit in 
a so-called continuous stranger sheet which is an orthotropic medium 
with finite extensional stiffness in the longitudinal direction, 
infinite extensional stiffness in the transverse direction, and finite 
shear stiffness. 

In order to investigate the possibility of a departure from the 
trend for a single cut, the analysis of reference 4 has been extended 
to two of the simplest double-cut cases. In each case the two cuts are 
separated by a single intact filament, with totals of two and four fila- 
ments being broken in the two cases. The dynamic results have been 
obtained by extending the single-cut dynamic analysis of reference 4- in 
the same way that the static analysis is extended. The Laplace trans- 
form is applied to equations (7) and (8), and initial conditions (9) are 
employed to obtain an expression for the transformed stress concentra- 
tion factor. Then a term-by-term inversion of the transformed stress 
concentration factor is carried out to obtain a series expression for 
the dynamic stress concentration factor which can be made as accurate 
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as desired by the retention of a sufficient number of terms . The 
dynamic analysis is presented in Appendix B. The reader is referred 
to reference 4 for a complete discussion of the analytical technique. 

The results are plotted in figure 7? where some of the results 
from reference 4 are reproduced for purposes of comparison. By 
comparison of the peak dynamic double-cut "values with the appropriate 
static -values from Table 1, it is found that the dynamic response 
factors for the cases of two and four broken filaments are 1.22 and 
1.23, respectively. This result is in keeping with the trend noted 
in reference 4 for single cuts. The fact that additional double-cut 
calculations have not been made because of their rapidly increasing 
complexity precludes the drawing of sweeping conclusions concerning 
the overall behavior of the double -cut dynamic response factor, 
nevertheless, radical departure from the single-cut trend appears to 
be unlikely. At the least, the present results do nothing to contradict 
the conclusion reached in reference 4 that dynamic effects appear to be 
of secondary importance in composite materials of the type investigated. 

Filament and Matrix Loads 

In this section, expressions are obtained for the most severe 
matrix shear forces for some single-cut and double-cut cases, in order 
to illustrate their behavior in the neighborhood of the cuts, and how 
they vary along the filaments. In addition, loads in broken filaments 
are calculated for some single-cut cases, and with the aid of an 
idealized model their load-recovery characteristics are examined in 
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Figure J. Dynamic stress concentration factors 
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connection with a current composite tensile failure model (see 
reference 8) . 

Matrix shear loads .- The dimensionless shear force per unit length 
between the n^ 1 and (n+l)^ filaments is given by 

s a U) = U n+1 (l) - u n (5) (37) 

In either a single-cut or double-cut configuration, the greatest matrix 
shear load is between a broken filament and an unbroken one. In the 
case of a double cut, the unbroken filament is an interior one. To 
find the matrix shear forces in a single-cut problem, the solutions to 
equations (24) for a specified number of broken filaments are substituted 
into equation (23), yielding the displacements of the filaments. Then 
the appropriate displacements are substituted into equation (37) to 
obtain the matrix shear load as a function of the longitudinal distance 
from the cut. For a double cut either equations (A4) or (A10) are 
solved, depending on the number of interior filaments, and their 
solutions are substituted into the second of equations (A3) or (A9), 
respectively, to obtain displacements which, m turn, are substituted 
into equation (37) to obtain the matrix shear load of interest. The 
influence functions V n ( i) required in equations (23), (A3) , or (A9) 
are given by equation (l6) , 

The results for single cuts across one, two, and three filaments 


are: 



38 


Single cut across one filament : 


S ° (S) = f 71 + 2 (^ (2S) ‘ Z o^ - |C L -1 (2|> - I l (2S) ) 


(58) 


where S n (0) = - 


° v ' 4 

Single cut across two filaments : 


s i<«> - f 


i^ 0 (ai) - MaDj + + “||l_ 2_(2!) - x 1 (£5^ - 


12 


(39) 


where S 1 (o) = - 

o 

Single cut across three filaments ; 


s 2 (£) = - 


3* 
1 6 


222)( M 2i) . ll(26) ] - ± - 


- L 0 ( 2 ^ 

+ i (20 

£ V 

4 - 74 

600 

3 t| Jt |3 

’ ' 


237 

£ 2 


m 


where S o (0) = - — rt. 

2 32 

The results for double cuts are: 

Two cuts, each across one filament, with, one interior filament : 


8 0 ( S )=i| 


-Is * - ^o(26)| + |fl + - Ix(2|)| 


_4_ _ JL2 


where S o (0) = ^ Jt. 

56 


m 
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Two cuts, each across two filaments, with one interior filament ; 


s 0 (l) 


5565* 

2(6512) 


22 + 



i 0 <2£5 - 





L_ x (2g) 



_4_ 


12 

*6? 


' I 4 1 + 7§ + S 


+ |)t-l (2£) - 1 1 (2£) ) - 7i ' 


;i 0 (24) 

36 _ 240 ^ 

jc|5y 


, 0 (2 £ )| ♦ |^3 + if 


(42) 


^ 5565 

where Sfo) = ~~.~-. z - - c-.. jt and the hars are used to denote double-cut 

0 2 ( 6512 ) 

shear loads. A comparison of equations (4l) and (42) will indicate 
how rapidly the complexity of double -cut shear force calculations 
increases with cut length. 

These most severe shear loads are plotted xn normalized form in 
figure 8 Also shown in figure 8 for easy reference is a table of 
their maximum amplitudes. As can he seen, the normalized shear loads 
are distinguished from one another mainly by their rates of decay with 
distance along the filaments, the decay being more gradual for longer 
cuts. It might also he noted that the peak shear loads for double cuts 
are less severe than those for single cuts across comparable or even 
smaller numbers of filaments. This would appear to be due to the fact 
that double cuts are directly adjacent to larger numbers of intact 
filaments (three or four) than are single cuts (always two)} hence, 
the load lost by the broken filaments is distributed over a larger 
number of intact filaments m the case of double cuts. 
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Figure 8. Variation of most severe matrix shear forces with axial distance from the cut 


Loads in broken filaments ." The calculation of loads xn “broken 
filaments xs of xnterest because it facilxtates an examination of the 
ability of the composite matrix to transfer the load back into the 
broken filaments and, hence, gives a measure of the degree of local- 
ization of the perturbed stress state. 

For a single cut across n consecutive filaments (starting with 
the zeroeth filament ) , the load in the i**- filament is given by 


n-1 


PiU) = 1 + Y %_ m ( 1)^(0) 

/ i 


m=0 


where the broken filaments are identified by i =0,1, •«•, n - 1. Each 

uu(o) is one of a unique set associated with a specific value of n- 

♦ 

The influence function W •( |) is given by equation (17) • Loads in the 
broken filaments have been calculated for a single cut across one, two 
and three filaments. The results are: 

One broken filament : 

P 0 (i) = 1 - - ^(21^ (^3) 

where P D (o) = 0. 

Two broken filament s : 


P D U) = pjU) = 1 - 






m 


where P o (0) = Pq(0) = 0 



Three broken filaments: 



Pi(i) = i 




Ix(2g) - L-i(2g)> 


+ 


2 L o (2S) - I 0 (2g) 




where P Q (0) = P^O) = P 2 (o) = 0. 

The results are plotted m figure 9* -She 90 percent load-recovery 
level is noted on figure 9 for later reference. It can he seen that 
the longitudinal distance required by a broken 'filament to recover a 
given fraction of its far-field load increases considerably with cut t 
length. This result is of interest in connection with the statistical 
failure analysis of reference 8, m which a significant parameter is 
the so-called "ineffective length" of a broken filament. The ineffective 
length is defined there as that portion of the length of a broken 
filament over which it supports less than 90 percent of its share of 
the load. In the next section, an idealized model is employed for 
calculating ineffective lengths with the aid of figure 9 > and. the 
implications of the results for the statistical failure analysis of 
reference 8 are briefly discussed. 
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Figure 10. Variation of ineffective filament length near a single cut 


Ineffective length calculation .- Studies cited in references 6 
and 8 suggest that the ineffective length, based on a load- recovery 
fraction of 90 percent of a single broken filament embedded in a matrix, 
can vary from one up to several hundred filament diameters, depending 
on the geometry of the composite and the mechanical properties of its 
constituents. Calculations have been made in reference 8, for example, 
of ineffective length as a function of filament volume fraction (ratio 
of filament volume to composite volume), Vf, and E/G, the ratio of 
filament Young's modulus to matrix shear modulus. The results are 
presented in the form of a family of curves, each member of which 
corresponds to a particular volume fraction. 

In the case of a composite containing a cut , an additional parameter 

which would be expected to influence the ineffective length is the 

number of broken filaments. With the aid of an idealized model, the 

present results can be used to obtain an indication of the influence 

of this parameter. A typical cross section of this model is shown in 

the sketch. For the present calculations, it is assumed that the 

thickness of the sheet is equal to the filament diameter. Then for 

jrdf 

this model the filament volume fraction is given by v^ = — where 
df is the filament diameter and d is the filament spacing. 

In terms of the pertinent mechanical and geometrical parameters, 
the axial distance from the cut is given by 
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x _ it E 


(V7) 


In order to obtain estimates of ineffective length changes due to the 

breaking of additional filaments, it is necessary to determine the 

values of £, denoted by £*, at which the broken filaments of interest 

have recovered 90 percent of their far-field load. These values of 

along with the appropriate values of E/G and v f , are then substituted 

into equation (47) to obtain values of the nondimen sional ineffective 

length, denoted by — . 

d f 
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This procedure has been carried out for two values of and a 

range of values of e/g. The cases of cuts across one, two and three 
filaments have been considered In the case of a cut across three 
filaments, the 90 percent recovery figure applies to either of the 
outer filaments m the cut rather than the middle filament, which would 
have yielded even greater values of ineffective length. The results 
are presented in figure 10 

The mathematical model employed in reference 8 consists of a single 
filament encased in a thin layer of shear-carrying material (binder) 
which, m turn, is embedded m an infinite body to which are assigned 
the average stiffness properties of the composite material. Hence, 
interaction between neighboring filaments is ignored. In addition, 
shear stresses in the average material are assumed to decay in a 
negligible distance from its interface with the thin layer of binder 
material. Because of these basic differences, the results of 
reference 8 are not included m figure 10. 

It should be noted, however, that the results of reference 8 yield 
a single curve for each filament volume fraction, whereas the present 
results yield a family of curves, each member of which corresponds to 
a different number of broken filaments. As can be seen, large changes 
in ineffective filament length can result from varying the number of 
broken filaments The failure analysis of reference 8, however, depends 
on a fixed value of ineffective length, regardless of the number and 
distribution of breaks m the composites. The present results suggest 
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that changes in ineffective length due to flaw size should he accounted 
for m statistical studies of the ultimate strength of filamentary 

composites. It might also he noted that the present results predict 

' l 

a more gradual recovery of load hy broken filaments than does the 
analysis of reference 8. 



CONCLUDING- REMARKS 


Analyses have been conducted of the loads and deformations in a 

i 1 

filament- stiffened sheet weakened by two collinear cuts and by periodic 

collinear cuts.' Also., some results have been obtained for the case of 

a single cut in addition to those previously reported. It has been 
, .*• 

found that significant interaction between collinear cuts is largely 
restricted to cases in which the distance between cuts is no greater 

j * * * 

than the cut length. It is seen that two closely spaced cuts can cause 
greater stress concentration than a single cut across a comparable total 
number of filaments. In the case of periodic collinear cuts., exact 
results indicate that a stress concentration factor derived for a 
transverse cut in an elastic strip in an earlier approximate analysis 
by Dixon (reference 3) is adequate for practical applications, although 
it is somewhat unconservative for cut lengths which approach the strip 
width. 

Limited calculations of dynamic stress concentration factors for 
suddenly introduced collinear cuts support the earlier conclusion that 
dynamic effects are not of great importance in filamentary composites 
of the type investigated. 

Maximum shear force calculations for single cuts of various 
lengths, which show that maximum shear forces grow more rapidly with 
cut length than maximum tensile forces, might be useful in determining 
whether a damaged composite material is more susceptible to shear 


failure or tensile failure. 
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Calculations of loads in broken filaments show significant changes 
m ineffective filament length with cut length, indicating that 
statistical strength analyses of composite materials should consider 
the incorporation of a flaw- size parameter. 

The present analyses are based on linear, small-deflection theory 
of a single filamentary sheet, whereas filamentary composites usually 

r 

are many filaments thick, and are subject to various nonlinear effects, 

t ** 

such as plastic deformations, large deflections and straightening of 
the filaments. Therefore, it appears that future analytical studies 
might be fruitfully directed toward analysis of better models of 
composite materials. 

* l 
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APPENDIX A 


EVALUATION OF INFLUENCE FUNCTION V n (|) 


The influence function V n (£) is given by 


* 1 r* 

V n (‘S) = ~ / cos 

Jiv 0 


n e e" 2 ^ sin 9 / 2 d0 


(Al) 


which can be changed to 


V (|) = - f cos 2n9 e“ 2 ^ sin 9 
n 3IW 0 


d0 


Let | 0; then (A2) becomes 


v = I f 11 (e 2in9 + e“ 2in0 )e 2:L $ sin 9 d0 

2itJo 


(A2) 


or 


- 1-J ^cos(2® sm0+ 2n0) + cos(2f sm 0 - 2n0) doj 


nit r 

- / £ 
2nJ 0 __ 


sin(2$ sin 0 + 2n0) + sin(2$ sin 0 - 2n0) d0 


0 
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which, "becomes 



J-2n(2*) * J2n(2^) 


[" 


sin(2n6 + 2$ sm 0) - sin(2n6 - 2® sin 0) 46 (A3) 


- 

) ■ 


where is the Bessel function of the first kind. The remaining 

integral may he evaluated to give 





J_2n ( 2 ®) + J 2n( 2 ®) 


+ ||B2n(-2») - E 2n (2$ j 


(A4) 


where E^ is the Weber function (see reference l) . Wow 
E 2ri (-2$)= -E„2n( 2 ^) 5 also, E_ 2n (2$) = E 2n (2$) for integer values 
of n, and J-2n( 2 $) = J ?n (2$) . Then equation (a 4) becomes 

V n = J2n( 2 ®) “ ^2n( 2$ ) ( A 5) 

Noting that $ = £e ijt / 2 , it is seen that 

JgnCa#) = (-l) n i 2n (2i) (A6) 

where I 2n is the modified Bessel function of the first kind From 


reference 1, for n an integer. 
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r 


~ H 0 (2$). 


n - 0 


E 2n< 2 ®> = \ n-X r f k + lW2n-2k-l 

$l \\ ' -W»>. 




k=0 r(2n + | - *) 


n = 1 


where H p^ is the Struve function. Making the change of variable gives 


r 


E 2n (2|el5t/2) = < 


- iL 0 (2£)> n = 0 


i(-D n+1 




L 2n (2S) + 


i y ^-l) » 
^ , 1 


2n-2k-l 


k=0 r| 


L_ 


(2n + | - kj 


-j 


, n > 1 


(A7) 


where L2 n is the modified Struve function (see reference l) . 
Substitution of (a 6) and (A7) into (A5) gives 

kr(k + | ) |2n-2k-l 

l2n(20 - W 2 ^) - ^ Z, ( “ 1) " 


V n (5) = (~l) n 


k=0 


'(an + | - fc) 


(A8) 


where the finite sum is to he omitted for n - 0. 
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APEENDIX B 

STRESS CONCENTRATION FACTORS FOR TWO 
EQUAL-IENGTH COELINEAR CUTS 


la this appendix attention Is restricted to problems in which the 
two collinear cuts are of equal length. This restriction is effected 
by setting q equal to s in equations (ll) and (12), which yields 

"N 

“(m-HL) r+s 

P n (6)=l- £ * X H n-i< 6 > U i< 0 > 

i=-(m+s) 


-(m+1) 


i=r+l 


r+s 




!!„(;) - i + Y, V n-ih> D i(°) + Y, Vi(i> D l(0) 

i=-(m+s) 


i=r+l 


J 


(Bl) 


and 


-(m+1) 


r+s 


° k =i+ Y; * Vi(°)^(o) - Y 

i=-(m+s) 


I=r+1 


(m + s) s n = - (m + l), r + l = n<r + s (B2) 


One of two cases arises, depending on whether the number of intact 
filaments between the cuts (interior filaments) is odd or even. The 

+ i. 

case of an odd number of interior filaments is treated first. 

Odd number of interior filaments .- When the two cuts are of equal 
length, the displacements are symmetric with respect to the line which 
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equally divides the group of interior filaments. In this case it is 
convenient to take the zeroth filament as the line of symmetry, so that 
m = r in equations (Bl) and (B2) , and the total number of interior 
filaments is 2r + 1. With 1^(0) = IL^O) because of symmetry, 
equations (Bl) and (B2) become after some manipulation 

i 

r+s * 

E n U) ' 1+ X /w n+ 1 (g) + N nil (i) j ujo)! 
i=r+l V ' 

r+s 

D nU) * ' + X 

i=r+l 


r 


n+i 


id) + Vid) u i(°) 




and 


r+s 

0 = 1+ Y k+i ( °) + k-i( 0) } U i^ 0 ^ r + l<n<r + s 
i=r+l ^ ^ 


(Bk) 


The maximum load will occur adjacent to the cuts (l = 0) m the two 
outermost interior filaments (n = lr). Hence, from the first of 
equations (B3) , the stress concentration factor for two cuts across 
s filaments separated by 2r + 1 filaments is 

r+s 

K s,2r-a = M°> = 1 + X 

i=r+l 

where the U^(0) * s are calculated from equation (B4) for specified 
values of r and s. 


r+i 


i(0) + Hr-i(O) %(0) 


(B5) 
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Even number of interior filament s . - In -fetus case, the i nterior 
filaments are assumed to be identified by - m ^ n < r, which with 
r = m - 1 becomes -m = n= m- l. Equations (Bl) and (B2) take the 
form 


-(m+l) 


m+s-1 


? n (0 =1+ ^ N n _iU)Ui(0) + y Vi(5) u i(°) 

i=-(m+s) ' • ‘ i=m 


-(m+l) ' m+s-1 

U n (i) = 5 + V Vn-iCDu^o) + Y v^OOu^o) 
, ' 5 ^-(m+s) “ 


) (B6) 


i=mi 


J 


and 


- (m+l) 


m+s-1 


°‘= 1+ . T . 

i=-(m+s) ' > i=m 


(m + s) s n = - (m + l), m = nSm+s-l (B7) 


< „ < 


For this configuration, the axial line of symmetry is located between 
the filaments numbered 0 and -1, which means that 

U -i(°) = u i-l(°) W 

After substitution from equation (b 8) and some manipulation, 
equations (B6) and (B7) become 
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r+s-l 


i=r 


r+s-l 


{* 




P n (i) = i + > K+1+1U) +Vi(£) K (o) 




(B9) 


u n (i) = i + 


. * O— JL. 

If 

i=r 


V n+i+l^ + Vx<*> K (0) 




and 


0 = 1 + 


r+s-l _ 

X l N n+i+l(°) + w n-x(°)}ux(0). 


r5n5r+s-l 


(BIO) 


In this case the maximum load occurs for n = - r, r - 1. Then with 

i 

n = r - 1 in the first of equations (B9 ) , the stress concentration 
factor for two cuts across s filaments/ separated ty 2r filaments, 
is 

r+s-l 

K s,2r = Pr-l(O) = 1 + X 

i-r 

where the U-^Oj’s are calculated from equations (BIO) for specified 
values of r and s. 


Hr+iCo) + N r _ 1 _ 1 (0)b- L (0) 


(Bll) 



APPENDIX C 


GENERAL PROOF OP FORMULA (2 6 )* 


From the system of linear algebraic equations 
n-1 

£ %_ r (0)U r (0) = - 1, 0 = i = n - 1 

r=0 

it is desired to obtain an expression for Uq(o) vhich is valid for 
all positive values of n. It is convenient to let 


W r = - u r (-o) 

x at 



and 


m “ n - 1 

Then the system (ci) becomes 
m 

X H i-k (0)W k = - *■> 0 = 1 = m 

. k=0 


^Communicated to the author by W. J. Harrington, Professor of 
Mathematics, North Carolina State University. 
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or, in an equivalent form 


r 


(-I-I)Wq + (l - + . 

( x - j) W 0 + + • 


• + 


v 2m - 1 2m + l y 


1 l > 

'W m = - 2 


. . + ( — I i — -)w = - 2 

V 2m - 3 2m - 1 m 


} (C2) 


i 2m - 3 2m - 


W 0 + * 




, 2m - 1 “ 2m + lj 0 


Wr> + » . . + 


• + (-l-DVl + f 1 ' jW - - 2 
- fW-1 + (-l-l)W m - - 2__ 


The system (C2) can he replaced hy a new system composed of the 
equations of (C2) linearly combined m the following manner: 

(1) First equation multiplied by (- ^ ) y 

(2) Sum of first two equations multiplied by 

(3) Sum of first three equations multiplied by (- i) j 

\ 

(m+1). Sum of first (m+1) equations multiplied by (- — -) . How let 

¥ k ^ + 2 
Wfc = — - i , and the new system has the form 
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The object is to calculate wq* To this end it would be advantageous 
to have a row vector which 'is orthogonal to all column vectors of the 
matrix A except the first- On page 102 of reference 7 is given the 
identity 



where 8 jjj 0 is the Kronecker delta. With (ie.,n = 2), this 

provides a row vector, * 
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which is orthogonal to the second column of the matrix A. It can he 
shown that this vector is orthogonal to all the columns of A except 
the first. The proof is by induction By way of example, it will he 
shown that the identity (C4) implies that 


m 

« Y 

k=0 




1 

2k - 3 


= 0 


if m ^ 2 


(C 6) 


Equation (c6) can he written 



(C7) 


Also, 


/2k + 2\ /2k \ / 2k \ 

y k + 1 j = 2 y k y + 2 ^k + lj 



1 

k + X 


Thus 

/2k 4- 2\ /2k\ / 2k\ 

I 2k + l) . 4 1 k ) 2l kJ 

2k - 1 (2k -*l)(k + 1) 



_ J L + k — L. 

3 2k - 1 3 k + 1 


2k - 1 
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Therefore, (C7) becomes 


& 5 (m) = - 



m-1 


■f | Gl(m - 1) + | 



'an - 2k - 2)\ 
m - k - 1 


k + 1 


(08) 


k.0 


On page 121 of reference 7, problem ll(a), appears an identity such 
that 



(C9) 


Thus if m > 1, then Gj_(m - l) =0, and from (c8) and (C9) there 
results 



Similar results can he obtained for G^(m) , G^Cm) , etc., from which the 
pattern of the induction proof can he recognized. 

Thus, for each m ? 1, the row vector (C 5) is orthogonal to all 
column vectors of the matrix A except the first, premultxplx cation of 
(C3) hy the row vector (C5) yields 



f 2k\(2m. 




(CIO) 


In problem 10 (a), page 121, of reference 7, the summation on the left- 
hand side of (CIO) is evaluated.; the result is 
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2 4m 


(2m + l) 



2 4 m (m .«) 2 

(2m + l).‘ 


Also, on page 130 of reference 7 the right-hand side of (CIO) is 
evaluated as 


(Cll) 



k=0 




(C12) 


By use of (Cll) and (C12), the solution to equation (CIO) is found to "be 


w 


o 



(2m * 1) 1 
2 2m (m l) 2 


or, replacing m by n - 1, 



(2n - 1).* 


2 2n-2Q n - i)ij 


n = 1,2,3, 


•which is equivalent to equation (26). 
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APPENDIX D 


SUMMATION OF THE INFINITE SERIES ^ W Rk+j(°) 

k=-co 

* * i ] 

The infinite series to "be summed is given by 


A(R,o) 




k=-co 


(DI) 


where %k+j(0) is given by 


W Rk+j(°) “ ~7~. T2 a 

Jt{4(KSt + j)^ - l} 


and where R is a positive integer, k is an integer, and j is an 
integer such that R > j . Then equation (Dl) can be written as 


A(R,j) 


kV 1 

4(Rk + j) 2 - 1 

k=-°° 


which can be changed to 


A(R,j) 



(D2) 


Therefore the problem is equivalent to summing the series 



66 


,,(a ’ rt ■ Z — a)V - 1) 


(D3) 


where a M L 0 [j + — )> a and 0 not integers. , 

R\ 2/ R\ 2/ 

The problem is amenable to a method presented m reference 2. We 


investigate the contour integral 



it cot tfZ 
(z - a)(z - 0) 


dz 


(D4) 


where Cq is the square contour with corners at z = + 

where m = 1, 2, 3> and k, and q is a positive integer- The presence 
of cot Jtz m the integrand of (b 4) provides a pole at each integer 
value of z. On each side of the contour cot it z is hounded, so that 
Inn Bq = 0. Calculation of residues at the poles enclosed by the 

q-»» 

contour then yields 


l\ ijt(2m.+l)/4 


lim Bq 
q-»» 


I cot jta cot 
jtl + 


a - 0 


r !)•! 

' k=- co 


(k - a)(k - 0) 


= 0 


/ i = — - — (cot itB - cot ita) 

L ' (k - a)(k - 0) a - 0 ' 


(D5) 


so that 
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Substitution of (C5) into (C2) , noting that a 
3 = - |(j + yields 



and 


X(R,j) 




or 
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APPENDIX E 

CALCULATION OF DYNAMIC STRESS CONCENTRATION FACTORS 


In dimensionless form the governing equations are 
Conservation of momentum: 

8 2 U n ct 2 U n . . 

— | + U n+1 - 2U- + U n _! = —£ * (ED 

d| 2 St 2 

r * • 

Boundary conditions: 

, i 

U (0,r) =0, - m = n = r, - (m + q + l) = n, n=x+s+l 

n 

P n (0,T) =0, - (m + q) = n = - (m + l), r + l = n= r + s (E2) 

P = 1 (E5) 

Initial conditions: 


P n U,o) = i 

8u n 

-^(6,0) = 0 


(e4) 


Application of the Laplace transform to equations (El) to (E3 ) and use 
of equations (E4) yield 


~n2 * 

Su n 

Bi 2 


+ <«. - (a + « 2 K + U S-1 - - » 


CB5) 
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with 


u *(0'£) ~ - m = n = r, - (m + q) ^ n, n^r+s+1 

= °> - (m + q) = n ~ - (m + l), r+l^n=r+s (E 6) 

- i 


where the first of equations (E^) has been converted to the condition 
U n U,0) = £, and the asterisks denote transformed quantities. 

Following the procedure of reference 4, use can be made of the 
unit solution to write the transformed loads and displacements m the 
form 

-(m+1) 

Pj(S,£)=i+ Y, N a-i(S>£K<°«£> 

5 i=-(m+q) 

r+s 

+ X N n-i^^) u i(°^) (E7) 

i=r+l 

-(m+l) 

u^U^)=|+ Yl V n-i(S,S)U*(0,5) 

i=-(m+q) 

r+s 

+ I Vj_i(5^)uj(0^) (E8) 

i=r+l 
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with the remaining boundary conditions (E6) talcing the form 


-(m+l) 

°=^+ Y, 

b i=-(m+q) 


r+s 

+ Y N n-i(°^)U^(0^>, 

i=r+l 


- (m + q) = n = - (m + l), 
r + ltsn^r + s 


(E9) 


Here the investigation is restricted to the cases of two equal- 
length collmear cuts separated by a single intact filament, which for 
convenience is taken to be the zeroeth filament. Thus, with m = r =0, 
q = s, and symmetry of displacements accounted for, equations (®7) and 
(E9) become 

Pj(t,E) =| + X foa(t*b + <-i(S>?)}t£(0,£) ' (E10) 

i=l 

and 


- i *Y (a 


-i < < 

1 = n = s 


(Ell) 


so that the cuts are separated by the zeroeth filament, and each cut 
severs s filaments. The transformed load in the zeroeth filament, 
which is the most highly stressed, is given by 

Pj(«,£) «| + Y + N* x U,£)|uf(0,£) (E12) 
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The stress concentration factor is the maximum value of this load, which 
occurs at £ = 0. Hence, the transformed stress concentration factor 
for two cuts across n filaments and separated hy a single intact 
filament is given hy 



i(C) 


*S(o,5) 


n 



!%(0,£) + Hli(0,£)>Ui(0,£) 


(E13) 


The transformed influence function was found in reference 4 to "be 





cos m0 



sin^ — + d0 
2 


(El4) 


From equation (El4) it can he seen that H m (0, £) = N_ m (0,£), so that 
equations (E13) and (Ell) can he written as 

n 

kJ,i(0 *= j + 2 £ nJ( 0,£)U*(o,£) ( El 5) 

1=1 

and 



For the cases in which each cut is across one and two filaments, 

* 7 

solution of equations (El6) for n = 1 and n = 2 yields 



• ff A 


(E17) 
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and 



%(N 0 + %) + N 2 (ff 0 + H 2 ) - (% + N 2 )(II 1 + N 5 )^ 

ug + 4 )(h; + <) - of + ^ 


(ei8) 


where , for brevity , the functional dependence of the transformed 
influence functions has been omitted. 

The inversion integral for the stress concentration factor is 


oS+loo «, 

Vi (T) - sr L „ 


d£ 


6-lco 


(E19) 


A series evaluation of this integral can be made by use of the method 
employed m reference 4 for the single-cut problem. Briefly, the 
method employs a conformal mapping function given by £ = z - l/z to 
transform the inversion integral (E19) into an integral around a 
contour C just inside the unit circle m the z -plane , given by 


t-iW -a/„ 


r 


a£,i 


z^ + l (z-1/z)t dz 

— e' — 

z 2 - 1 z 


(E20) 


This integral is evaluated by finding the coefficient of the zeroeth 
power of z in the expansion of 

_ Jz d - 1 


( E21) 
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The -berm in brackets in (E21) can be expanded m a power series m 
by using equation (E17) or (El8) in conjunction with the series 
expansion for Nj^(0,z), which was found in reference ^ to be 

CO 


5^(0, z) = - 


1 

z 



/ Q ,.2k+m 

(-z2) 


z 


(E22) 




and 


k — 1,2,3, * • 


■X - 

Since the expansion for £K n ^i involves only positive even powers of 
z, only the negative even powers m the expansion of the other part 
of (E21) need be sought. It was found m reference k that if 

prrf e <*‘i/ a)T -Z V* 

k=-eo 

then 


Cq — - 1 

c_ k = J 0 (2 T) + 2J 2 (2t) + . . . 2J k _ 2 (2r) + J k (2r) - 1, (E23) 


k = 2, 4, 6, . 
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after considerable manipulation,, the results are found to be 

k i,i ( t) " 1 - c a (T) + k c s ( t) + i c io< T > - \ °l^ T> 

2 

+ C l8( T ) + C 22( t ) + ‘ • 

2 y 2 


and 


K2 A (t) = 1 - C 2 (t) - \ 0 4 (t) - | C 6 (t) + 1 C 8 (t) + C 10 (t) 


15 , . 49 . , 8i . . 47 . . 

+ ^ T) + p? C ^ (T) + pi c i 6 (t) + pi c is‘ T > 


- C 20 (t) - 2^2 C 22 (r) + . . 


>15 


,17 


(E24) 


Good convergence is obtained for the range of r investigated. 



